Magnetoelastic coupling in RETiQ 3 (RE = La, Nd, Sm, Gd, Y) 
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A detailed analysis of the crystal structure in i?£Ti03 with RE — La, Nd, Sm, Gd, and Y 
reveals an intrinsic coupling between orbital degrees of freedom and the lattice which cannot be 
fully attributed to the structural deformation arising from bond-length mismatch. The TiOe oc- 
tahedra in this series are all irregular with the shape of the distortion depending on the RE ionic 
radius. These octahedron distortions vary more strongly with temperature than the tilt and ro- 
tation angles. Around the Ti magnetic ordering all compounds exhibit strong anomalies in the 
thermal-expansion coefficients, these anomalies exhibit opposite signs for the antiferromagnetic and 
ferromagnetic compounds. Furthermore the strongest effects are observed in the materials close to 
the magnetic cross-over from antiferromagnetic to ferromagnetic order. 
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I. INTRODUCTION 



The rare-earth titanates (RET1O3) with a distorted 
perovskite structure have a single electron in the t2 S or- 
bitals of the Ti-3d shell and are all Mott insulators^. 
Therefore, they can be considered as the one-electron 
counterpart to the cuprates with a single hole in the 
3d shell. Upon decreasing the ionic radius in the rare- 
earth RE series from La to Y— , the calculated one- 
electron bandwidth gets smaller— and the deviation of 
the Ti-O-Ti bond angle from 180° increases 4 -. Concomi- 
tantly, the Ti magnetic order changes from G-type an- 
tiferromagnetism with the magnetic moment in a di- 
rection (RE = La,...,Sm) to ferromagnetism with the 
magnetic moment in c direction (RE=Gd,...,Yb and 
Y) 5 i 6 i 7 i 8 iV°i u i 12 i 13 i 14 ' 15 i 16 . These apparently distinct 
magnetic ordering schemes, however, belong to the same 
irreducible representation, as it has been clarified long 
ago^i for the structure type of GdFeC>3 to which all 
i?_ETi03 belong. Due to the low symmetry of the crys- 
tal structure, the G-type antiferromagnetic component 
in a direction is always coupled with an A-type antifer- 
romagnetic moment in b direction and a ferromagnetic 
component in c direction. Thus the change from anti- 
ferromagnetic to ferromagnetic order of the Ti moments 
in the series of RETiOz is just a redistribution of the 
ordered moment between the three different components 
(G x , Ay and F z ) within the same magnetic symmetry. 

Tremendous efforts have been made to understand 
the magnetic properties of these titanates during the 
last years^^ 1 ^^! 2 ^ 2 ^. An orbital liquid model 
has been proposed in the idea that the splitting be- 
tween the t-ig states would be sufficiently low to allow 
orbital fluctuations to play a dominant role in the phys- 
ical propertie d 20 ' 21 . On the other hand, in a more con- 
ventional scenario, it is argued that orbital fluctuations 
would be efficiently suppressed by the t 2ff -level splitting 
of the order of 200 meV^l^^L. Taking account of 
the anisotropic orbital arrangement Schmitz and Miiller- 
Hartmann could explain nearly all the details of the mag- 
netic order in LaTiOs including the measured dispersion 



of the magnetic excitation s 24 ' 25 . 

In addition to the Ti magnetism, magnetic order of 
the RE moments occurs at lower temperature. For the 
lighter RE ions (RE = Ce, ... , Sm) the RE moments 
order antiferromagnetically below Tff2 < Tni (with Tjvi 
being the Neel temperature for G^-type antiferromag- 
netic ordering at the Ti sites). For the i?i?Ti03 com- 
pounds studied in this work belonging to this first group 
(RE = Nd, Sm) an antiferromagnetic ordering of F y /C z - 
type is observed at the RE sites§i For the RET1O3 with 
heavier RE ions (RE = Gd,...,Yb) there are two basic 
magnetic structures. In the first class, the RE magnetic 
moments order ferromagnetically and couple antiparallel 
to the ferromagnetic component of the Ti sites (F z ). In 
the second class (RE=Tb, Dy and Ho) the magnetic mo- 
ments at the RE sites have both ferromagnetic and an- 
tiferromagnetic components in the a, b pian o 16 ' 28 ' 29 ' 30 ' 31 . 
Among the heavier RET1O3 materials only GdTi03 has 
been studied in this work, it exhibits ferromagnetic Gd- 
moment ordering (first group). These complex RE mag- 
netic structures, however, seem to yield little insight to 
the Ti magnetism and the associated interaction ener- 
gies are smaller than those in-between the Ti's except in 
the compounds where the Ti-Ti interaction is strongly 
suppressed. 

The GdFe03-type structure of the rare-earth titanates 
differs from the ideal, cubic perovskite structure (with 
space group Pm3m) by a rotation of the octahedra 
around the c axis and by an additional tilt of the TiO"6 
octahedra around the cubic [110] axis. While the oc- 
tahedral tilts are opposite for all neighboring TiOg oc- 
tahedra, the octahedra are rotated in the same sense 
along the c direction and opposite only in the ab planes. 
The combination of tilt and rotation distortion reduces 
the symmetry to the orthorhombic space group Pbnm 
and enlarges the unit cell to a \/2 x y/2 x 2 cell ro- 
tated by 45° around the c axis (see Fig[T]). In this space 
group internal distortions of the TiOg octahedra can oc- 
cur without a further symmetry breaking. The most fa- 
mous example for such additional distortions is LaMn03 
where the antiferroorbital order of the e g electrons co- 
exists with the tilt and the rotatio n 32 ' 33 ' 34 . Qualita- 
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FIG. 1: (color online) Crystal structure of RET1O3 with 
corner sharing Ti06 octahedra in Pbnm notation. Dashed 
lines show the cubic unit cell of the undistorted structure and 
solid lines show the orthorhombic unit cell. Green: Ti-ions, 
blue: O-ions and red: RE-ions. 



tively similar octahedron distortions have been observed 
for all rare earth titanates studied in this work^ 5 . For 
an undistorted octahedron, one expects a negative or- 
thorhombic splitting s = (a — b)/(a + b) since the octa- 
hedral tilt takes place around the orthorhombic b axis 
(Pbnm). But due to an elongation of the TiOg octahe- 
draon basal plane along the orthorhombic a direction the 
orthorhombic splitting in LaTi03 becomes even positive. 
In the series of i?£Ti03 the type of octahedral distortion 
changes its character with decreasing rare-earth ionic ra- 
dius. For LaTiOa the 02-02 distances (01 denotes the 
apical oxygen displaced from the Ti sites along c and 02 
the basal in-plane oxygen) along the octahedron edges 
differ strongest resulting in a rectangular shape of the 
basal oxygen planes. In contrast, in YTi03 there is lit- 
tle difference in these 02-02 distances but the two in- 
plane Ti-02 distances split. These distortions together 
with the anisotropic Ti-RE coordination lead to a split- 
ting of the t2 g energy levels with a different character of 
the lowest state. Interestingly the orbital arrangement 



resulting from these two distortions is totally different. 
Within a single layer parallel to the a, b plane, the varia- 
tion of the edge lengths yields a ferroorbital configuration 
whereas the splitting of the Ti-02 distances yields an an- 
tiferroorbital configuration 5 . According to Goodenough- 
Kanamori rules, these distinct orbital ordering schemes 
agree with the antiferromagnetic and ferromagnetic Ti 
ordering in LaTi03 and in YTi03, respectively. 

Although there is no doubt that the octahedron dis- 
tortions in the RET1O3 series are coupled with the mag- 
netic ground state, it is less obvious whether these effects 
should be considered as resulting from orbital ordering. 
There is no symmetry breaking in the titanate case as the 
orbital distortions are embedded in the crystal structure 
deformation arising from tilt and rotation. Pavarini et 
alj22 argue that the values of the tilting and the concomi- 
tant shift of the RE sites would determine the splitting 
of the orbital levels^. However, the lattice parameters 
and the octahedron deformation in LaTi03 exhibit clear 
anomalies around the magnetic ordering suggesting that 
there is an intrinsic coupling between orbital degrees of 
freedom and the crystal structure 5 -^. It is the aim of 
the present work to further elucidate this orbital-lattice 
coupling in the RET1O3 series and to separate the intrin- 
sically electronic effects from tilt and rotation distortions. 
For this purpose we have analyzed the known crystal 
structure of oxides showing the same GdFe03 structure 
type but having no electron in the 3d shell. Furthermore 
we have studied the temperature dependence of the crys- 
tal structure in i?_ETi03 for RE = La, Sm, Nd, Gd, and 
Y in great detail. In all compounds we find evidence for a 
direct orbital-lattice coupling in addition to that arising 
from the distortions. 



II. EXPERIMENTAL 

Single crystals of RET1O3 were grown using a floating- 
zone image furnace, as described elsewhere^. Parts of 
them were crushed in order to obtain powder samples of 
high quality. Magnetic ordering temperatures of the Ti 
ions were determined by a SQUID or a vibrating sample 
magnetometer. The magnetic ordering of the Ti depends 
sensitively on the stoichiometry REi^^TiOs+s- In order 
to obtain crystals with a high transition temperature it 
is necessary to tune the parameters during the growth 
process. In this work several crystals of the same mate- 
rial were grown and analysed; the highest Tjv/Tc-values, 
we find, agree with or are higher than those reported in 
a recent study 16 . The Neel temperatures of the anti- 
ferromagnetic compounds are 146 K for LaTiOs; 94 K 
and 81 K for two samples of NdTi0 3 and 54 K and 48 
K for two crystals of SmTi03. In the case of NdTi03 
the high T^r-sample was used for diffraction and ther- 
mal expansion studies and the lower Tjv-sample only for 
diffraction studies. For SmTi03, the 54 K sample was 
studied in thermal expansion and the 48 K sample by 
diffraction techniques. The Curie temperatures for the 
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ferromagnetic titanates are 36 K (RE=Gd) and 29 K 
and 25 K for two YTi03 samples used for diffraction and 
thermal expansion measurements, respectively. The Sev- 
eral diffraction techniques were applied to analyze the 
temperature dependencies of the structural parameters. 
A Siemens D5000 diffractometer equipped with a one- 
dimensional detector was used for powder X-ray diffrac- 
tion with Cu-K Q or Cr-K Q radiation. These as well as 
the neutron powder diffraction data were analyzed by 
the Rietveld technique, but the laboratory X-ray pow- 
der diffraction does not allow one to obtain the oxygen 
structural parameters. In order to perform full struc- 
ture analyzes powder neutron-diffraction patterns were 
recorded on the high-resolution 3T.2 diffractometer at 
the Orphee reactor in Saclay (A = 1.225lA). The full 
crystal structure is also determined by single-crystal X- 
ray diffraction studies on a Bruker X8 Apex CCD diffrac- 
tometer using Mo-K Q radiation. Absorbtion correction 
for spherical samples was performed with the standard 
Bruker software (scale) and isotropic extinction cor- 
rection was applied during the refinement on F 2 with 
Jcma20 00 38 . Some of the crystals were twinned. Due 
to very similar lattice parameter in LaTi03, the compo- 
nents of all twin domains could be integrated together 
and refined with different volume fractions for all twin 
domains within Jana2000. The corresponding twin laws 
result from a threefold rotation about the [111] axis of 
the 8 A pseudocube combined with a reflection across the 
(1-10) plane. In the other i?i?Ti0 3 , no twinning was ob- 
servable or the contributions of the different twins could 
be separated during integration procedure thus obtaining 
a dataset only for the major twin domain. This was the 
case in NdTi03, where the reflections of different twin do- 
mains which could not be separated from each other were 
dismissed (see also Tab. Q}. These minor twin domains 
have only small volume fractions of roughly 10% of the 
total sample volume. The structural results of our sin- 
gle crystal X-ray and powder neutron-diffraction experi- 
ments are summarized in Tab. I & II. Thermal-expansion 
coefficients were measured for single-crystalline YTi03, 
SmTi03, and NdTi0 3 on a high-resolution capacitance 
dilatometer— . For YTi0 3 we furthermore measured the 
magnetostriction along the a direction in magnetic fields 
up to 14 T. 



III. RESULTS AND DISCUSSION 

A. Analysis of known GdFeOs-type structures 

Since it is not obvious in the RET1O3 series to separate 
the intrinsic orbital-ordering effects from those imposed 
by the severe tilt and rotation distortions, it appeared 
interesting to analyze the published crystal structures 
of metal oxides AMO3 in GdFe03-type structure which 
posses an empty 3c? shell. Using the Inorganic Crystal 
Structure Database (ICSD)^ we have studied the com- 
pounds : CaTi0 3 , CaZr0 3 , EuSc0 3 , GdGa0 3 , GdSc0 3 , 
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FIG. 2: (color online) (a) The octahedral tilt O and the ro- 
tation $ versus the tolerance factor t for various d° systems 
with perovskite structure (space group number 62, Pbnrn^). 
(b) Difference of both basal M-02 bond lengths and difference 
of both 02-02 bond lengths versus t. (Open symbols denote 
the values for LaTiOa and YTKD3; lines are linear fits to the 
data.) 



HoA10 3 , LaGa0 3 , NdGa0 3 , PrGa0 3 , SmA10 3 , SrZr0 3 , 
TbA10 3 with respect to their tilting (0) and rotation ($) 
distortions and to the deformation of the metal-oxygen 
octahedron. 4 - We discuss the relative difference between 
the two M-02 distances, A(M — 02), and that between 
the two 02-02 basal-plane edge lengths, A((92 - 02). 
The analysis of the published structures is resumed in 
Fig. O The tilt and rotation deformations arise from the 
bond-length mismatch between the M-0 and A-0 bonds 
which in the case of an ideal perovskite would be forced 
into a ratio of y/2. With the ionic radii 2 one calculates 
the tolerance factor t = (r^ + ro)/[V2~ • (tm + r o)] with 
values much below t = 1 indicating a strong mismatch 
due to too small A ionic radii. Fig. [2] (a) shows the 
tilt and rotation angles as a function of the tolerance 
factor yielding the expected increase of both distortions 
towards smaller tolerance factors. Also the ratio between 
the tilt and rotation angles seems to be fully determined 
by the tolerance factor. One may note that LaTi0 3 and 
YTi0 3 match these relations. In Fig. [2] (b) we show the 
A(M — 02) and A((32 — 02) values again as a function of 
the tolerance factor. For the A(M — 02) bond-distance 
ratio there is no clear tolerance-factor dependence visi- 
ble, but the deformation of about 2.5% seen in YTi0 3 is 
much larger than those in all the others materials which 
vary only between 0% and 1%. This gives the first evi- 
dence that the A(M — 02) deformation is not purely a 
consequence of the tilt and the rotation in the titanates. 
Concerning the edge-length deformation characterized by 
A(02 — 02), the deformation of about 3% seen in LaTi0 3 
is strongest but there seems to be a sizeable distortion 
also in the series without 3d electrons. This deforma- 
tion is intimately coupled with the orthorhombic split- 
ting. Tilting of a rigid octahedron around b will always 
yield a smaller a lattice parameter, whereas some of the 
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materials, including LaTiC>3, show the opposite. This re- 
versed orthorhombic splitting, however, appears only for 
small tilt angles. It seems to result from the interactions 
between the A site and the apical oxygens. One of these 
bond distances is significantly shortened by a rigid tilt. 
The systems may therefore elongate perpendicular to the 
tilt axis in order to reduce this effect. This will yield a 
weak elongation of the octahedron basal plane along a 
and eventually the reversed orthorhombic splitting. Al- 
though LaTiC>3 exhibits the strongest A(02 — 02) devi- 
ation, additional information is needed in order to clarify 
an underlying orbital effect. 



B. Temperature dependence of the crystal 
structure in antiferromagnetic RETiOs (RE = La, 
Nd, Sm ) 

For the RE titanatcs with larger rare-earth ionic radius 
showing antiferromagnetic order, we observe a strong rise 
of the a lattice parameter at the onset of the magnetic 
ordering on cooling. Concomitantly, the b lattice param- 
eter decreases more rapidly leading to a strong increase 
in the orthorhombic splitting e = (a — b)/(a + b) upon 
cooling. This anomalous behaviour of the a lattice pa- 
rameter was first observed for LaTiO^^. Our measure- 
ments on other RET1O3 show that this anomalous rise of 
e becomes even enhanced for smaller RE in the series of 
i?_ETi03 as long as the system stays antiferromagnetic, 
see Fig. [He). 

In Fig. [4]the results of powder neutron diffraction mea- 
surements on LaTiOs up to T=750 K are shown. The 
structure data at the first three temperatures were taken 
from Ref £ where the same sample as in the present study 
was used. Over the full temperature range, the variation 
of both the octahedral tilt O and the rotation $ is very 
small, but the ratio of the 02-02 octahedron edge lengths 
changes essentially, from approximately 4% at 10 K to 
approximately half of this value at 750 K. Thus only the 
internal octahedral distortions vary strongly with tem- 
perature. The decoupling of tilt and rotation angles on 
the one side and of the octahedron distortion on the other 
side demonstrates that the latter is not just a conse- 
quence of the former. The structural bond- length mis- 
match does not drive the octahedral distortion. Instead, 
the structural anomaly has to be ascribed to an orbital 
effect which, however, does not break any symmetry but 
just enhances the crystal-field splitting promoted through 
the tilt and rotation distortions. 

Since the octahedral tilt O is nearly temperature in- 
dependent, the rise of the orthorhombic splitting e can 
be taken as an indicator for the elongation of the Ti06 
basal plane A (02 — 02), because only O and this distor- 
tion define the value of e. Indeed, the plot of A (02 — 02) 
versus e shows a linear dependence (see Fig. 

As can be seen in Fig. [3J the low-temperature increase 
of e is more pronounced for the i?i?Ti03 with smaller 
RE, NdTi03 and SmTi03, still showing antiferromag- 
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FIG. 3: (color online) (a-c) Orthorhombic lattice parameters 
of RETi0 3 (RE=L&, Nd, Sm, Gd and Y) as a function of 
temperature, (d) c parameter scaled to the room temperature 
value, (e) Orthorhombic splitting e — (a — b)/(a + b). The 
arrows indicate the magnetic transition temperatures Tjvi / 
To and the stars mark the values for e measured by powder 
neutron diffraction. Open circles denote the values for the 
second Ndi-^TiOa+a sample with a different stoichiometry 
and Tjv = 81 K. 



netic order. Furthermore, with decreasing rare-earth 
ionic radius the rise of e starts closer to the magnetic 
ordering temperature, even though it sets in well above 
Tn for all samples studied. In SmTi03 with a Tn of only 
45 K the anomalies are most pronounced. The increase 
of e by 0.01 between 100 K and 10 K indicates an en- 
hancement of the octahedral distortion characterized by 
A (02 — 02) by 1% , which is quite remarkable for an 
internal structure parameter. For SmTiOs the tempera- 
ture dependence of the crystal structure has been studied 
by single-crystal X-ray diffraction down to 100 K. Again 
there is almost no variation in the rotation and tilt an- 
gles whereas the internal octahedral deformation signif- 
icantly changes. Compared to LaTi03 the increase of 
A(02 - 02) between 300 and 100 K is much smaller but 
most likely much stronger between 100 K and lowest tern- 
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FIG. 4: (color online) Results of our powder neutron diffrac- 
tion measurements on LaTiOs at different temperatures. The 
first three points were taken from Ref.— . The others were 
measured on the same sample, (a) TiC>6 octahedral tilt O 
and rotation $. (b) Differences of the 02-02 bond lengths, 
(c) Differences of the Ti-O bond lengths, (d) Differences of 
the 02-02 bond lengths versus e. The approximately linear 
dependence of A(02-02 bond length) and e underlines that 
£ is a good indicator for the octahedral distortion A(02-02 
bond length). (Lines are guide to the eyes.) 



peratures, see Table HI Qualitatively, the observed struc- 
tural anomalies at the antiferromagnetic ordering also 
agree with a mechanism based on spin-orbit coupling 42 , 
since the elongation of the octahedron occurs parallel to 
the ordered spin moment pointing along the a direction, 
but the observed effect appears to be too large. 

As can be seen in Fig. 3, we have studied a second sam- 
ple of Ndi-^TiOa+a, which due to minor inhomogeneity 
exhibits a lower Neel transition at XV=81 K. Clearly, 
the structural anomalies are strongly suppressed in this 
material. 

In LaTiOs no structural anomaly is seen along the c 
direction, see also Ref<2£, but NdTiC-3 and SmTiC>3 ex- 
hibit such anomalies at the Neel ordering. Upon passing 
into the antiferromagnetic state the lattice expands along 
the c direction. The relative strength of these anomalies 
is, however, still much lower than those appearing in the 
ab planes. Since the ground-state orbital 5 -* 2 ^ in LaTiC>3 
has significant extension along c this elongation agrees 
still with the orbital-ordering effect. 



C. Temperature dependence of the crystal 
structure in Ferri-/Ferromagnetic RET1O3 (RE = 
Gd, Y) 

The temperature dependences of the lattice constants 
a and b in GdTiC>3 and in YTi03 are totally different 
from that in the titanates with antiferromagnetic or- 
dering discussed above. The ferromagnetic compounds 
show a negative thermal-expansion coefficient along the 
b direction and a larger positive one along a. Again, 
there are well-defined anomalies around the magnetic or- 
dering, which, however, exhibit the opposite signs com- 
pared to the antiferromagnetic compounds. The anoma- 
lies are visible in Fig. [3] In the case of YTiC>3 these 
low-temperature structural anomalies were already re- 
ported in reference^ 5 -. The anomalies are, however, much 
stronger for GdTi03. 

The crystal structure of YTi03 was studied by neu- 
tron diffraction on a powder sample at 35 K, i.e. slightly 
above Tq, and at 2 K and by single-crystal X-ray diffrac- 
tion on a different sample between 100 K and room tem- 
perature. Besides the changes in the lattice constants 
discussed above, no significant structural changes could 
be detected between 100 K and room temperature. Fur- 
thermore, the powder data at 35 and 2 K did not reveal a 
measurable change in the crystal structure as well. The 
Ti-02 bond-length distortion is rather large in YTi03, 
A(M-02)=2.91% (at 100 K), whereas the other param- 
eter is small, A(02 - O2)=1.08% (at 100 K). Again, the 
rotation and the tilt angles vary very little between room 
temperature and 2 K. The structural anomalies observed 
in these ferromagnetic titanates cannot be explained by 
the spin-orbit coupling. The most prominent octahe- 
dron elongation occurs along the b direction, whereas one 
would expect the spin-orbit coupling to cause an elonga- 
tion along the c direction which is the direction of the 
ferromagnetic ordered moment 42 . Note, however, that 
indeed there is such an elongation along the c direction 
upon passing in the ferromagnetic state. In addition to 
the spin-orbit coupling, the orbital-ordering effect may 
also contribute to the c-axis expansion, as the lowest or- 
bital level in YTi03 also possesses significant extension 
along o±&±&22. 



D. Magnetostriction and thermal-expansion 
measurements 



In order to better resolve the structural anomalies, we 
have measured the thermal-expansion coefficients a = 
1/L ■ dL/dT for three RETi0 3 compounds with RE = 
Nd, Sm and Y (see Fig. [5|) on a high-resolution capac- 
itance dilatometer— . The results are given in Fig. [5] 
and should be compared to the similar measurement on 
LaTiO-3 described in reference 35 . For these studies we 
have used almost fully untwinned single crystals, as it can 
be ascertained by the comparison between the dilatome- 
ter and the diffraction results. 
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As expected from the diffraction results, the thermal- 
expansion anomalies in SmTi03 and in NdTi03 are much 
stronger than those in LaTiC-3; in particular the anomaly 
in the c parameter of NdTiOs is well resolved, whereas 
it has no counterpart in the high-resolution studies on 
LaTi03 35 . The high-resolution studies in SmTi03 and 
in NdTiOs confirm our conclusion that the structural 
anomalies set in far above the Neel ordering^. Further- 
more, the smaller anomalies at the ferromagnetic order- 
ing in YTiC>3 are confirmed by the dilatometer study. 
In YTiC>3, the c-axis anomaly is larger and the anomaly 
along the a direction is reversed in comparison to the 
antifcrromagnetic compounds. 

In Fig. [6] we show the results of magnetostriction ex- 
periments on YTi03. The measurement was performed 
approximately along the a-directionii The anomaly in 
the thermal-expansion coefficient broadens and shifts to 
higher temperature with increasing magnetic field. This 
field dependence is expected for a ferromagnetic transi- 
tion. Applying the magnetic field at low temperature well 
in the ferromagnetic phase, we first observe an elonga- 
tion along the a direction of about — ~ 1-10 -5 . At zero 
field the magnetic moments in YTiC>3 order ferromagnet- 
ically along the c axis; under an external field along the 
a direction they first flip parallel to the field. The small 
elongation at small field can hence be associated with the 
reorientation of the ordered moment. This elongation is 
most likely caused by the spin-orbit coupling elongat- 
ing the octahedron along the ordered moment^. This 
low-field magnetostriction is, however, smaller than the 
magnetostriction occurring at higher field. The high-field 
effects have to be considered as an enhancement of the 
integrated zero-field thermal-expansion anomalies, as it 
can be seen in Fig. [5] a). At zero field the ferromagnetic 
ordering in YTiC>3 is accompanied by a small antiferro- 
magnetic moment along the a axis due to the competing 
magnetic interactions^, which are determined through 
the orbital arrangement. Upon cooling in zero field, the 
orbital arrangement is changing in the way that ferro- 
magnetism is further stabilized; this effect seems to be 
further strengthened by the external field which anyway 
stabilizes ferromagnetic order. 



E. Crystal structure across the RET1O3 series 

Fig. [7] shows the magnetic phase diagram of i?_ETi03 
as a function of the rare-earth ionic radius. In addition 
to the pure RET1O3 samples we have included results on 
samples where the A site is occupied by a mixture of La 
and Y. The diagram clearly shows how the crossover from 
ferromagnetic to antiferromagnetic order is driven by the 
rare-earth ionic radius, even though some influence of ad- 
ditional effects is visible. In the ferromagnetic RET1O3, 
there is a stronger variation of the Curie temperature 
which indicates the direct influence of the RE-Ti mag- 
netic interaction. Non-magnetic Y implies a comparably 
low Curie temperature in this series. Nevertheless, the 
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FIG. 5: (color online) Thermal expansion a. = 1/1- dl/dT for 
(a) NdTiOs parallel to b and c, (b) SmTiOs parallel to a and 
(c) YTiC>3 parallel to a and c. 



rare-earth ionic radius can be considered as the main ex- 
ternal parameter driving the magnetic transition in the 
pure as well as in the mixed compounds. 

In Fig. [8] the results of our single crystal diffraction 
measurements performed at room temperature are re- 
sumed. Qualitatively, these results agree with previous 
studies^ but there are significant quantitative differences 
due to the improved sample quality and due to he higher 
precision rendered possible with the modern CCD-X-ray 
techniques. 

Both NdTiC-3 crystals with T at of 94 K and of 81 K 
were studied at room temperature. The lower Tn of the 
second crystal already indicates some non-stoichiometry. 
With X-ray diffraction it is difficult to determine pre- 
cisely the oxygen content, therefore, it is not astonishing 
that the change in the refined oxygen content per for- 
mula unit from 2.952(30) (high-T^ sample) to 2.988(12) 
(low-T^v sample) is not very significant. However, in the 
low-Tjv sample we find a significant amount of vaccan- 
cies on the Nd-site with a content per formula unit of 
only 0.979(1), whereas the high-T^y sample as well as the 
other RET1O3 compounds with high transition temper- 
atures show the ideal RE to Ti ratio within comparable 
precision. The structural parameters obtained for the 
low-Tjv and high-Tjv samples are comparable in view 
of the large variation of the internal parameters in the 
i?_ETi03 series, see Table I, however clear differences can 
be detected. The main influence of the excess oxygens 
consists in a reduction of the octahedron tilting. Since 
excess oxygen and Nd- vacancies enhance the Ti valence 
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FIG. 6: (color online) (a) Thermal expansion a — l/l ■ dl/dT 
of YTiC>3 for various magnetic fields as a function of temper- 
ature (H\\a). (b) Magnetostriction AL(H)/L of YTi0 3 as a 
function of applied magnetic field H\\a. 



and thereby reduce the effective Ti-ionic radius, this tilt 
reduction is in perfect agreement with the bond-length 
mismatch scenario described by the tolerance factor. A 
similar suppression of the tilt distortion by excess oxy- 
gen was also found in La2Cu04+a 4 -. The enhanced Ti- 
valency in the non-stoichiometric sample is directly ob- 
served in the bond-valence sum which is 0.045(1) larger 
than that in the stoichiometric compound. 

In Fig. [5] we plot the structural parameters against 
the rare-earth ionic radius. The largest overall structural 
changes concern the tilt and the rotation deformations. 
The decrease of the rare-earth ionic radius enhances the 
bond-length mismatch and these angles, as seen in Fig. 
M(a) and in the Ti-O-Ti angles shown in Fig. [£](&). The 
distortion of the ideal perovskite structure gets extremely 
strong for YTiC>3 which appears close to the stability of 
the Pbnm structure. Due to increasing tilt and rotation 
deformations the RE-0 distances vary strongly within 
the series; the ratio of the longest to the shortest bond 
increases from 1.42 in LaTiC>3 to 1.63 in YTiC>3. For 
a rigid tilt an 01 site directly moves towards a RE site 
yielding the shortest RE-0 bond in the entire series (RE- 
O distance #6 in Fig. [8]) . This short bond seems to cause 
the inversed orthorhombic splitting seen in many weakly 
distorted perovskites, as discussed above. For the smaller 
rare earths, however, the systems seems to avoid the ex- 
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FIG. 7: (color online) Magnetic phase diagram for RET1O3 
and Lai-j; Ya;Ti03. Tc and Tjv are plotted as a function 
of the rare-earth ionic radius. Squares: Tc/Tjv for RETiOs 
from this work, triangles: values for .RETiOs from Greedan^, 
crosses: values for i?_BTi03 from Katsufuji et al.—, circles: 
values for Lai-^Y^TiOs from Okimoto et aZ.— , stars: values 
for Lai_^Ya;Ti03 from Goral et al^. The inserts show the 
principal octahedral deformations of the end members. 



tremely short RE-Ol bond distance by shifting the RE 
site. As can be seen in Fig. [8] (g) the shift of the RE site 
with respect to its high-symmetry position at (0.5,0,0.25) 
varies stronger than the oxygen internal parameters re- 
lated with the tilting. Due to this RE shift, the shortest 
RE-Ol distance does not get very small and at the same 
time the distance to the opposed 01 sites increases less, 
see Fig. |8] (c). Furthermore, the RE-site shift in y di- 
rection results in a strong splitting of the two distances 
to the two 01 sites near (0.5, ±0.5, 0.25). The structural 
change implied through the smaller ionic radius in the 
i?_ETi03 series is, therefore, not only described by the 
increasing tilt and rotation angles, but in addition there 
is a change in the RE-0 coordination. 

With the increasing shift of the RE site for smaller 
rare-earth ionic radius the deformation of the RE-Ti co- 
ordination increases as well, as it is shown in Fig. [S] 
(h). In the cubic perovskite, each RE has eight near- 
est Ti neighbors; these distances significantly split for 
the smaller RET1O3 most likely destabilizing the orbital 
arrangement of the LaTi03 type^^. 

For the Ti-orbital physics the internal parameters of 
the TiOe octahedra are important. In agreement with 
the previous discussion 5 there is a crossover in the octa- 
hedral distortion through the RET1O3 series. Whereas 
LaTi03 exhibits the rectangular elongation of the basal 
plane characterized by A (02 — 02), YTi03 exhibits 
the rhombic distortion characterized through distinct in- 
plane Ti-0 bond distances and A(M — 02). As stressed 
in reference^, this structural crossover is accompanied 
by the magnetic crossover from the antiferromagnetic to 
the ferromagnetic ground state. A part of this crossover 
seems to be induced by the structural effects discussed 
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FIG. 8: (color online) Results from single crystal X-ray 
diffraction measurements at room temperature (open circles), 
(a) Apical tilt and basal rotation $ of the TiOe octahedra 
of -R£Ti0 3 . (b) Ti-O-Ti bond angles, (c) RE-0 bond lengths 
(see plot (g) for notation), (d) 02-02 bond lengths, (e) Ti-O 
bond lengths, (f) ratio of 02-02 bond lengths and ratio of 
Ti-02 bond lengths (see also plot (h)). (g) Displacement of 
the RE ion from the symmetry position 00 i. (h) RE-Ti bond 
lengths, (i) Oxygen coordination of the RE i+ ion with the 
numbering for plot (c). (j) TiOe - black arrows indicate both 
02-02 bonds which are about parallel to a or b and white 
arrows indicate the two Ti-02 bonds which are rotated by $ 
towards a or —b. 



above, i.e. the response of the RE-0 arrangement to the 
extremely strong tilt and rotation distortions. For the 
smaller rare earths the displacement of the RE site is 
more important and has a stronger direct influence on the 
crystal-field splitting of the t2 g levels. This purely struc- 
tural explanation of part of the deformation is further 
supported by the similar effects occurring in the REFeO^ 



series, whose structure was studied in detai l 48 ' 49 i 5[l . The 
REFeO^ series shows almost the identical tilt and rota- 
tion angles as a function of rare-earth ionic radius causing 
similar constraints in the RE-0 coordination. Interest- 
ingly also this system exhibits the octahedral basal-plane 
elongation for larger rare earths and the Fe-02 bond dis- 
tance splitting for smaller rare earths, but both effects 
are much weaker than in the case of the corresponding 
titanates. 

With the analysis of the intrinsic structural response 
to the smaller rare-earth ionic radius we may more deeply 
discuss the structural anomalies occurring in all i?_ETiC>3 
studied as a function of temperature. It appears that the 
tilt, the rotation and the RE-site shifts implied through 
the rare-earth ionic radius already favor a certain crystal- 
field splitting 22 ' 23 , which in case of LaTiC>3 results in 
a ferroorbital ordering (within the afr-plane) whereas it 
causes an antiferroorbital ordering in YTiC>3. The com- 
pounds close to the structural crossover, where the nature 
of the lowest orbital level should change as well, appear 
most interesting. The ferroorbital and the antiferroor- 
bital ordering imply an antiferromagnetic and a ferro- 
magnetic nearest-neighbor magnetic interaction, respec- 
tively. Close to the structural and magnetic crossover 
the magnetic interaction is hence significantly reduced 
in agreement with the lower values of Tq and Tn (see 
Fig. [7]). Upon cooling, each of the titanates passes into a 
well ordered state and in order to suppress the magnetic 
fluctuations the magnetic interaction parameters should 
get enhanced on the cost of a lattice deformation. Since 
either the ferromagnetic or the antiferromagnetic interac- 
tion should increase upon cooling, either the antiferroor- 
bital (RE= Gd, Y) or the ferroorbital (RE— La, Nd, Sm) 
ordering (within the a6-plane) must get enhanced by ro- 
tating the orbital degree of freedom. This explains first 
the sign change of the in-plane anomalies in the series 
and second the observation of the strongest anomalies 
in the compounds close to the crossover. Upon continu- 
ously varying the magnetic nearest-neighbor interaction 
J, it has to pass zero which due to the diverging magnetic 
fluctuations yields an energetically unfavorable phase. In 
all compounds, the orbital effects follow the crystal-field 
splitting implied by the structural deformation, thereby 
enhancing the magnetic interaction. 

As suggested in the phase diagram shown in Fig. [7] the 
rare-earth ionic radius should be considered as the exter- 
nal parameter driving the magnetic transition from anti- 
ferromagnetic to ferromagnetic ordering. Therefore, the 
i?_ETi03 phase diagram may be compared to the many 
known systems where a magnetic transition may be im- 
plied by varying some external control parameter. As it 
has been recently discussed for a generic quantum phase 
transition one may expect a sign change in the thermal- 
expansion anomalies when crossing the transition^ 3 - due 
to the fact that the system will avoid the accumulation 
of entropy just at the transition. The sign changes in the 
thermal-expansion anomalies were for example observed 
across the metamagnetic transition in Ca2- a; Sr a ;Ru04 54 . 
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F. La^Y.TiOs (x=0.50) 

Since the rare-earth ionic radius drives the structural 
deformations, the orbital arrangement and the mag- 
netic structure, it would be quite interesting to study 
this phase diagram by continuously varying the ionic 
radius. La^-rY^TiC^ was thought to be a good sys- 
tem since there is no magnetic moment on the RE site. 
Several groups have already studied this system, see 
reference s 51 ' 52 . Indeed this system shows the transition 
from antiferromagnetism to ferromagnetism between 30 
and 40 % of Y— . We have prepared a powder sample 
of Lao.5oYo.5oTi03 and show the inverse of its magnetic 
susceptibility in Fig. [5] (b). There is no well-defined 
anomaly comparable to the effects in the pure i?_ETiC>3 
but the weak effect around 90 K indicates the onset of 
antiferromagnetic ordering in qualitativ agreement with 
the Curie- Weiss analysis of the susceptibility at higher 
temperatures. From these magnetic properties and from 
the averaged ionic radius which is between those of Nd 
and Sm one would expect similar physical properties as 
those observed in NdTi0 3 and in SmTi0 3 . Fig. M (a) 
shows, however, that this is not the case. Although the 
orthorhombic lattice parameters are in between the corre- 
sponding values of the pure Nd and Sm compounds, com- 
parable structural anomalies at low temperatures do not 
occur in Lao.soYo.soTiOs, see the plot of the orthorhom- 
bic distortion in Fig. [5] (a) . We think that the disorder 
arising from the mixed occupation of La and Y with quite 
different ionic radii destroys the effects induced by orbital 
ordering. Although this result annihilates the hope to 
tune a mixed titanate to the critical averaged rare-earth 
ionic radius and to drive it thereby into a critical configu- 
ration, it corroborates our conclusion that the structural 
anomalies in the RET1O3 materials arise from orbital 
physics. 



IV. CONCLUSION. 

The detailed structural analysis of the RET1O3 series 
as a function of temperature allows us to separate the 
structural effects implied directly through the variation 
of the rare-earth ionic radius from those related with an 
orbital effect. The decrease of the rare-earth ionic radius 
leads first to the enhancement of the tilt and rotation 
angles and second to a modification of the RE-0 and RE- 
Ti coordinations. For small distortion angles the RE- 
site is closer to its high-symmetry position whereas it 
moves away from that for the larger distortions. As has 
been stated by several groups these shifts yield already 
a significant crystal field splitting of the t2 9 levels. The 
character of the thereby lowered orbital is rather different 
and the orbital arrangement, which is either ferroorbital 
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FIG. 9: (color online) (a) Orthorhombic splitting e of 
Lai-zY^TiOs (x=0.50) compared with the other i?£Ti0 3 
with similar rare earth ionic radii, (b) Inverse magnetic 
susceptibility 1/x of Lao.sYo.sTiOs compared with LaTiOa, 
YTiOs and the average of LaTiOs and YT1O3. Arrows indi- 
cate Tc/Tjv. 

or antiferroorbital in nature, is even opposed for the end- 
members of the series. 

The temperature dependence of the crystal structure, 
however, shows that anomalous structural deformations 
occur in all RET1O3 which should modify the char- 
acter of the orbital with the lowest energy. Since the 
ferroorbital and antiferroorbital a6-plane arrangement 
is related with an antiferromagnetic and ferromagnetic 
nearest-neighbor interaction, respectively, the structural 
change directly couples to the magnetism. The struc- 
tural anomalies are clearly associated with the magnetic 
ordering temperatures where the thermal-expansion 
coefficients exhibit extrema. We, therefore, conclude 
that magnetism drives a change in the orbital ground 
state. In all compounds, this change points in the 
sense to enhances the dominant magnetic interaction. 
For antiferromagnetic RET1O3 with La, Nd and Sm 
the ferroorbital arrangement is strengthened while for 
GdTiC>3 and YTiC>3 the antiferroorbital arrangement 
is strengthened. In this sense the temperature-driven 
orbital ordering follows the crystal-field splitting already 
imposed through the strong distortions. The structural 
effects are weakest in the end members and strongest 
close to the structural and magnetic crossover. Unfortu- 
nately, it seems not possible to drive the titanates into 
the critical configuration by continuously varying the 
ionic radius in a mixed system Lai-^Y^TiOs as disorder 
effects seem to suppress orbital degrees of freedom. 
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TABLE I: Results of single crystal X-ray diffraction measurements of RET1O3 (RE = La, Nd, Sm, Gd, Y) at room temperature 
and for RE — Sm, Y down to 100 K. *: results for a (non-stoichiometric) NdTiOa-sample with Tat = 81 K. }: The LaTiOa single 
crystal was twinned with the following ratios of the six possible domains: 82.7%, 8.1%, 4.2%, 1.8%, 1.5% and 1.5%. f: Both 
NdTiOa-crystals were slightly twinned (a <-» b) with a major twin domain of roughly 90% and all reflections with \h — k\ < 1 
(h, k 7^ 0) were excluded as they could not be separated from other twin contributions in the integration procedure. 
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R,E: 


La 


La 


Y 




17 


1 


550 


'4 ' 


2 
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M c ( MB /Ti 3 +) 








0.753(75) 
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20 


lattice: 












a (A) 


5.6386(2) 


5.6454(2) 


5.3226(1) 


5.3239(1) 
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b (A) 


5.6367(2) 


5.6447(2) 


5.6952(1) 


5.6944(1) 
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c (A) 


7.9518(2) 


7.9699(3) 


7.5962(2) 


7.5952(2) 




atoms: 
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x(REl) 


0.9928(5) 


0.99329(1) 


0.9776(2) 


0.9780(2) 




y(REl) 


0.0387(2) 


0.0353(3)) 


0.0740(1) 


0.0738(2) 




z(REl) 


0.25 


0.25 


0.25 


0.25 




x(Til) 
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y(Til) 


0.5 


0.5 


0.5 


0.5 




z(Til) 
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x(Ol) 


0.0778(5) 


0.0759(7) 


0.1213(2) 


0.1212(2) 


y(oi) 


0.4882(4) 


0.4885(5) 


0.4570(2) 


0.4572(2) 




z(Ol) 


0.25 


0.25 


0.25 


0.25 
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x(02) 


0.7102(4) 


0.7113(5) 


0.6901(1) 


0.6901(2) 


27 


y(Q2) 


0.2908(4) 


0.2896(5) 


0.3092(2) 


0.3095(2) 




z(02) 


0.0401(2) 


0.0391(3) 


0.0577(1) 


0.0576(1) 




B (A ): 












B(RE1) 


1.07(2) 


1.33(3) 


0.177(12) 


0.193(12) 




B(Til) 


0.55(3) 


0.70(4) 


0.190(25) 


0.181(26) 




B(Ol) 


1.05(3) 


1.32(4) 


0.257(16) 


0.305(17) 
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B(02) 


1.21(3) 


1.60(3) 


0.228(10) 


0.264(10) 


fit: 
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R 


4.31% 


4.45% 


3.61% 


3.70% 


R, 


5.78% 


5.93% 


4.89% 


4.91% 
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TABLE II: Results of the rietveld refinement of powder neu- 
tron diffraction measurements of RET1O3 (RE = La and Y) 
at different temperatures. 
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